Experimental evidence of chaotic advection in a convective flow 
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Lagrangian chaos is experimentally investigated in a convective flow by means of Particle Track- 
ing Velocimetry. The Finite Size Lyapunov Exponent analysis is applied to quantify dispersion 
properties at different scales. In the range of parameters of the experiment, Lagrangian motion is 
^ found to be chaotic. Moreover, the Lyapunov depends on the Rayleigh number as TZa 1 ^ 2 . A simple 

, dimensional argument for explaining the observed power law scaling is proposed. 

b ! ... , . „, 

The investigation of transport and mixing of passive tracers is of fundamental importance for many geophysical 
and engineering applications [j^] . It is now well established, and confirmed by several numerical j^] and experimental 
j| evidence, that even in very simple Eulerian flow (i.e. laminar flow) the motion of Lagrangian tracers can be very 
complex due to Lagrangian Chaos In such a situation, diffusion may be of little relevance for transport which 

is, on the contrary, strongly enhanced because of chaotic advection 

In this Letter, we address the problem of quantifying the dispersion of passive tracers in a relatively simple convective 
^| | flow at various Rayleigh numbers IZa. By applying the Finite Size Lyapunov Exponent (FSLE) Q analysis to the 
. Lagrangian trajectories obtained from Particle Tracking Velocimetry (PTV) technique, we are able to estimate the 
dispersion properties at different scales. We find a clear power law dependence of the Lagrangian Lyapunov exponent 
on the Rayleigh number. This dependence is explained by a dimensional argument which excludes a role of diffusion 
>— h ■ in the dispersion process. 

i " i l The experiment is performed in a rectangular tank L — 15.0cm wide, 10.4cm deep and H — 6.0cm height, 
filled with water. Upper and lower surfaces are kept at constant temperature while the side wall are adiabatic. 
The convection is generated by an electrical circular heather of radius 0.4 cm placed in the mid-line of the tank at 
0.4 cm above the lower surface. The heather works at constant heat flux Q which is controlled by a feedback on the 
power supply. By changing the heat input we control the Rayleigh number IZa = (g[3QH 3 ) /{olvk), where g is the 
gravitational acceleration, (3 the thermal expansion coefficient, a the thermal conductivity, v the kinematic viscosity 
and k the thermal diffusivity. In the parameters range explored in our experiments the flow consists of two main 
counter-rotating rolls divided by an ascending thermal plume above the heat source . The upper end of the plume 
t-H ■ oscillates horizontally almost periodically with a frequency which depends on the Rayleigh number. 

Lagrangian trajectories are obtained by PTV technique [^|. The fluid is seeded with a large number of small (50 /im 
in diameter), non- buoyant particles. The vertical plane in the middle of the tank and orthogonal to the heat source 
is illuminated by a thin laser light sheet. Single exposure images are taken by a CCD camera and subsequently 
digitalized at 8.33 Hz rate with a 752 x 576 pixels resolution. Trajectories are then identified as time ordered series 
of particle locations. 

Each run lasts for 2700 s, corresponding to 22500 frames. Typically 900 particles are simultaneously tracked for 
each frame. In the following we analyze the trajectories obtained in 6 different runs with Rayleigh number in the 
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range 6.87 10 7 < IZa < 2.17 10 9 . 

The analysis of Lagrangian data have been done with the Finite Size Lyapunov Exponent tool which, introduced 
in the context of the predictability problem turbulence |6j, has already been demonstrated very efficient for the 



characterization of dispersion in bounded domains [|9| and in the treatment of experimental data |10|. Let us recall 
the basic ideas on the FSLE; more details can be found in The idea of the Finite Size Lyapunov Exponent 

is to generalize the Lyapunov exponent, which measures the average rate of divergence of two infinitesimally close 
trajectories, to finite separations. To this aim, we fix a set of thresholds R n = Rop n (n = 0, . . . , N) and we consider, 
at each time t, new couples of trajectories Xi(i), X2(i) found at separation R(t) = |xi(£) — X2(i)| < Rq. We follow the 
evolution of the trajectories and compute the "doubling time" T p (R n ) it take for the separation to grow from scale 
R n up to R n +i = pR n - Of course it must be p > 1, but we take p not too large in order to avoid contributions from 
different scales. 

After performing a large number of doubling time experiments (over the possible different couples in the run) we 
average the doubling time at each scale R from which we define the Finite Size Lyapunov Exponent 

A(i?) = (r(k lnp ' « 



1 



where < [. . .] > e indicates the average performed on the doubling time experiments (see Refs. |6||| for further details). 
The FSLE is a generalization of the Lyapunov exponent A || , in the sense that 

lim X(R) = A , (2) 

fl->0 

physically speaking X(R) « A for R < Ie, where Ie is the smallest Eulerian characteristic length. For larger values of 
R, X{R) gives information on the mechanism governing the dispersion at scale R. For example, in the case of standard 
diffusion, on the basis of dimensional considerations one expects that 

\{R)~D/R 2 , (3) 

where D is diffusion coefficient. 

The use of the FSLE is particularly useful for studying the dispersion properties of passive tracers in closed basins 
and, therefore, fits very well with the considered flow. In such a situation, asymptotic regimes like diffusion (j^) 
might never been reached due to the presence of boundaries. Denoting by R ma x the average couple separation in the 
asymptotic uniform distribution, it has been found that for a large class of systems, for R close to Rmax, A(i?) fits 
the following universal behavior 

\{R)^^ Rma l R , (4) 

TR R 

where tr has the physical meaning of the characteristic time of relaxation to the uniform distribution. Equation (^) 
can be obtained assuming an exponential relaxation of tracers' concentration to the uniform distribution 

We applied the FSLE analysis to the Lagrangian trajectories experimentally obtained. In order to increases 
the statistics at large separations R, we have computed X(R) for different values of the smallest scale Ro (Ro = 
0.4cm, 0.6cm, 0.8cm). The threshold ratio is p = 1.2 for all the analysis. For the results presented below we use 
H = 6 cm as unit length and the diffusive time t K = H 2 /k ~ 25000 s as unit time. 

In Figure [l] we show the X(R) versus R computed for the run at TZa = 2.39 10 8 . The first important result is the 
convergence of X(R) to the constant value A ~ 3100 i" 1 at small R. This corresponds to an exponential divergence 
of close trajectories, i.e. a direct evidence of Lagrangian chaos. The large value of the Lyapunov exponent (in unit 
of inverse diffusive time) indicates that chaotic advection is the dominant mechanism for particle dispersion at small 
scales. 

For larger separation X(R) drops to smaller values, indicating a slowing down in the separation growth. This is 
quantitatively well described by the saturation regime ([!]). The collapse of the curves at different Rq confirms that 
sufficient high statistics is reached even at these large scales. Fluctuation among different Ro curves can be taken as 
an estimation of the error for X(R). By fitting the large scale behavior of A(i?) with (Q) we obtain R max ~ 1.9 H 
and tr ~ 8.010 _4 i K . Thus also for the late stage of relaxation to the uniform distribution, diffusion seems to play 
a marginal role. The characteristic Eulerian scale in the flow Ie can be estimated by the end of the exponential 
regime (plateau X(R) — A) at which the non linear effects start to dominate. We find Ie — 0.5 H, indeed not too far 
from the saturation value. In this condition there is no room for the development of a diffusive regime Q, as Fig. 1 
clearly shows. In addition, let us mention that by comparing X(R) computed at different TZa, we find that both the 
characteristic scales Ie and the saturation scale R ma x are independent on the Rayleigh number. 

In order to explore the dependence of the Lagrangian statistics on the Eulerian characteristics, we have performed 
the FSLE analysis for Rayleigh number varying over more than one order of magnitude. The dependence of the 
Lagrangian Lyapunov exponent A (computed from the plateau of A(i?) at small R) on TZa is shown in Figure |^. A 
clear scaling is observed, indicating a power law dependence 

A ~ TZa 1 (5) 

with 7 = 0.51 ± 0.02. 

It is worth noting that because of the geometry of our experiment, the flow shows qualitatively the same pattern 
for the whole range of TZa explored. This is confirmed by the independence of Ie and R ma x on IZa as discussed above. 
In these conditions, the scaling of A on TZa can be supported by the following dimensional argument. The equations 
of motion in the Boussinesq approximation, made non-dimensional by H and t K , are |llj|: 



where Pr = vJk is the Prandtl number (which is kept constant in our experiments). It is now easy to verify that 
performing the following rescaling: 



u a — ► Au a , t — > A t , TZa — ► A Ka , (8) 

where A is an arbitrary factor, equations (^0) remain unchanged, if the diffusive terms are neglected (as it is suggested 
by the previous results) . Let us stress that this rescaling do not involve neither the space (as it is suggested by above 
discussion) nor the Prandtl number (because in our experiments we change Q, keeping both v and k constant). The 
consequence of the Eulerian scaling invariance on the Lagrangian motion, governed by 

^ = u(x(t),t), (9) 

is that Lagrangian trajectories are independent on the Rayleigh number. The Lyapunov exponent, which is dimen- 
sionally the inverse of a time, thus scales as TZa 1 ' 2 , according to the result shown in Figure 2. 

As a consequence of the dimensional argument (§-§) the FSLE should scale with TZa both in the linear and in the 
saturation ranges, i.e. \{R) /TZa 1 / 2 is a i?a-independent function. In order to verify the scaling argument, in Figure 
H we plot the FSLE compensated with TZa 1 / 2 for the different runs. The collapse is fairly good, and the fluctuations 
are of the order of the error reported in Figure |2|. As a further argument in favor of the scaling (|^) , let us observe 
that the transitions to the saturation regime occurs at a scale independent on TZa. 

As a conclusion, let us stress again that the approximate agreement with the dimensional scaling argument and 
the large values of the Lyapunov exponent in units of diffusion time give a coherent picture of Lagrangian dispersion 
ruled by chaotic advection. 
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I. FIGURE CAPTIONS 
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FIG. 1. X(R) versus R for different initial thresholds R = 0.067// (o) , 0.1 if, (A), 0.13 H (y) at TZa = 2.39 10 8 . The 
straight line is the Lyapunov exponent A = 3100 ± 200 1^ 1 and the curve is the saturation regime (W) with tr — 8.0 10 -4 t K and 

Rmax — 1.9 if . 
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FIG 2. Lagrangian Lyapunov exponent dependence on the Rayleigh number TZa. The errors are estimated by the fluctua- 
tions corresponding to different initial Rq. The line is the best fit A ~ TZa ' 51 . 
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FIG. 3. Data collapse of X(R) at different IZa rescaled with TZa 1 ^ 2 . 
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